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Abstract 

By explicit construction of the ADHM data, we prove the existence of a charge 
seven instanton with icosahedral symmetry. By computing the holonomy of this 
instanton we obtain a Skyrme field which approximates the minimal energy charge 
seven Skyrmion. We also present a one parameter family of tetrahedrally symmetric 
instantons whose holonomy gives a family of Skyrme fields which models a Skyrmion 
scattering process, where seven well-separated Skyrmions collide to form the icosa- 
hedrally symmetric Skyrmion. 
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1 Introduction 



Skyrmions are a type of topological soliton in three-dimensional space which are 
of interest to physicists, in that they are candidates for a solitonic description of nuclei. 
Numerical simulations reveal that the minimal energy Skyrmions often have a great deal 
of symmetry, and in particular the charge seven Skyrmion has icosahedral symmetry [|J. 

One approach to the study of Skyrmions is the suggestion by Atiyah and MantonP] 
that low energy Skyrmions of charge n may be approximated by computing the holonomy 
of a charge n instanton in Euclidean M 4 along lines parallel to the Euclidean time-axis. 
This proposal therefore predicts the existence of a charge seven instanton with icosahedral 
symmetry. In this paper we verify that indeed such an instanton exists by presenting 
the corresponding ADHM data. We then use this data to compute the holonomy of the 
instanton and hence generate a Skyrme field of charge seven with icosahedral symmetry, 
which is an approximation to the minimal energy Skyrmion. 

Manton has proposed |16| that the low energy dynamics of n Skyrmions may be ap- 
proximated by motion on a finite dimensional manifold of charge n Skyrme fields, and one 
reasonable choice for such a manifold appears to be the moduli space of n-instantons[@]. 
In such an approximation a scattering of seven Skyrmions is described by a one-parameter 
family of 7-instantons. We present such a family of instantons, obtained by imposing 
tetrahedral symmetry, and compute the associated Skyrme fields. We display the baryon 
density isosurfaces for this scattering process, in which six Skyrmions approach along the 
Cartesian axes a seventh Skyrmion at the origin. The Skyrmions merge to form first a 
dodecahedron, then a cube, and finally the dual dodecahedron before again separating 
along the Cartesian axes leaving a single Skyrmion remaining at the origin. Part of this 
scattering process, where the dodecahedron deforms to its dual, is believed to be an im- 
portant vibrational mode of the dodecahedral Skyrmion, and these are of relevance when 
considering the quantization of Skyrmions||. 



2 An icosahedral 7-instanton 

In this section we shall present the ADHM data for an icosahedrally symmetric 7- 
instanton. First, we briefly recall the ADHM construction of instantons IJ and give an 
explanation of symmetry in a gauge theory. 

The ADHM data for an SU(2) n-instanton consists of a matrix 

Mm) <"> 

where L is a row of n quaternions and M is a symmetric n x n matrix of quaternions. By 
the ADHM constraints we shall mean the condition 

M f M is a real matrix, (2.2) 

where denotes the quaternionic conjugate transpose of M. 
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The first step in constructing the instanton from the ADHM data is to form the matrix 

A <*> =(*-,!.) < 2 ' 3 > 

where l n denotes the n x n identity matrix and x is the quaternion corresponding to a 
point in K 4 via x = x 4 + %X\ + jx 2 + kx 3 . 

The second step is then to find the (n + l)-component column vector N(x), of unit 
length, which solves the equation 

N(x)^A(x) = 0. (2.4) 

The final step is to compute the gauge potential A^(x) from N(x) using the relation 

A^x) = N(x)^d^N(x). (2.5) 

This defines a pure quaternion which can then be regarded as an element of su(2) using 
the standard representation of the quaternions in terms of the Pauli matrices. 

In order that these steps make sense it is necessary that the ADHM data satisfy an 
additional invertibility condition: that the columns of A(x) span an n-dimensional quater- 
nionic space, for every x. Equivalently this condition can be expressed as: 

A(x)^A(x) is invertible for every x. (2-6) 

In the absence of this condition ADHM data satisfying Q2.2|) only gives rise to a self-dual 
gauge field with singularities (corresponding to the points x where fl2.6|) fails). 



2.0.1 Symmetric instantons 

It will be useful to start by considering the problem of symmetric instantons in a general 
context. Thus let G C SO (3) be a subgroup of the rotation group of Euclidean 3-space, 
and make it act on M 4 by rotations on (xi, X2, £3), leaving x 4 alone. This action has a very 
convenient description using quaternions. For this, replace G by the corresponding binary 
group G C SU (2) (the double cover of G obtained from the double cover 577(2) — > 50(3)). 
Now 

SU(2) = Sp(l) = {u G H : uu* = u*u = 1} (2.7) 

and x 1— ► uxu~ l clearly preserves the £4 (real) component of x. It also acts on the imaginary 
part of x by the rotation in 50(3) corresponding to u in SU(2). (In fact this gives a 
construction of the double cover of 50(3) by SU(2).) From now on, whenever we speak 
of a subgroup of 50(3) acting on M 4 , we shall mean that it acts in the way we have just 
described. 

In order to get a grip on the problem of imposing symmetry in a gauge theory, let 
us introduce a rank-2 vector bundle E over R 4 and assume that our gauge potential 
defines a connection on E. If G acts on E then it makes good sense to require A^ to be 
invariant (or symmetric) under this action. In this case an action of G on E consists of 
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the data of a unitary map Q(g,x) : E x — > E gxg -i for each g G G and x G R , depending 
smoothly on x, and satisfying Q(l,x) = the identity, and 

Q(f, gxg-'Mg, x) = Q(fg, x) for all f,g,eG,xe R 4 (2.8) 

both sides of this equation being maps E x — > E^g^j^-i . 

It is important to note that one cannot assume that there exists a gauge in which 
all the f2's are equal to the identity. What is true, however, is the following. Let a be 
any fixed-point of the action of G, i.e. any point on the ^4-axis. For such a point, (|2.8| ) 



reads fi(/, a)£l(g } a) = £l(fg, a) for all /, g, G G, giving a representation of G on the fibre 
E a . This representation changes to an equivalent one under changes of gauge, so if this 
representation is non-trivial, there cannot be any gauge with the fi's all equal to the 
identity. 

Conversely, given any complex 2-dimensional representation p of G, we may define an 
action of G on a bundle over R 4 by putting 

Q(g,x)=p(g) (2.9) 

relative to some background gauge. Obviously this action restricts to the representation p 
on the fibre E a . In fact, if G is a finite subgroup of 577(2), any action of G on E is of this 
form for some representation p. (This is so because R 4 is topologically trivial and can be 
contracted to the fixed-point set of the action. In general the classification of actions on 
bundles is more complicated.) 

Given a bundle E with an action Q, we can state the condition for a gauge potential 
to be symmetric: 

A^(gxg- 1 ) = n(g, x)A^x)n(g, x)' 1 - d^g, x)n{g, x)~\ (2.10) 

In this context, the JTs are often referred to as 'compensating gauge transformations' a 
practice we shall occasionally follow in this paper. The term — <9 M f2(g, x)Q(g, x) _1 vanishes 
if the action is in the standard form ( |2.9| ). Notice that if fl2.10| ) is satisfied then the cur- 
vature (field-strength) transforms like F(gxg~ 1 ) = Q(g,x)F(x)Q(g,x)~ 1 and for example 
the action density —ti{F llv F ,J,v ) is a G-invariant function on R 
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Technical Remark Since the object of physical significance, the gauge potential, takes its 
values in the adjoint bundle of E, it is most natural to impose symmetries on this bundle. 
Thus we suppose that G acts on Ad(E), covering the action of G on H (so g G G gives a 
linear map Ad(E) x — > Ad(E) gxg -i for every x in H). Then it does not necessarily follow 
that such an action of G on Ad(E) lifts to an action of G on E. Moreover this failure 



occurs in at least one physically interesting example, the cubic instanton of 

To explain how this problem can be understood, we remark that, as above, the action of 
G on Ad(E) restricts to give a representation of G on Ad(E) a for any point a fixed by the 
action. Since Ad(E) a may be identified with the Lie algebra of SU{2) and hence with R 3 , 
in general the image of G by this representation will be a subgroup of SO (3) isomorphic 
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to some quotient group F = G/H of G. In general the double cover F C SU(2) of F acts 
on E a and this extends, as in Q2.9|) to an action on the whole of E. However it is in general 
not the case that F is a quotient of the binary group G. In the case of the cubic instanton, 
G is the rotation group of the cube, H is the Klein viergruppe (consisting of the identity 
and the half-turns about the three coordinate axes) and F is isomorphic to the dihedral 
group of order 6. 

In the general case, it still happens that a double cover G, say, acts on E, but this need 
not be the binary group G. With F, F as above, 

G = {(gJ)eGxF:p(g) = q(f)}. 

Here p : G — > G/H = F is the map to the factor group and q : F — > F is the double-cover. 
The restriction to G of the projection G x F — * G is then a 2 : 1 map, as one can easily 
verify from the definitions. 

If G is the icosahedral group, this complication does not arise owing to the fact that G 
is then a simple group (being isomorphic to the alternating group A§) and H must either 
be 1 or G. 



2.1 Symmetric ADHM data 



Returning to the ADHM description (2J3), this will be G-symmetric if for every g G G, we 
have 'compensating gauge transformations' 



A^- 1 ) = ^ j A(x)U(g)-' (2.11) 

for every <? in G. The first matrix has been decomposed into blocks corresponding to ( |2.1| ) 
so poo is 1 x 1 and U is n x n. In order to preserve the shape of A(x), U(g) must be the 
product of a real orthogonal matrix with a unit quaternion, while Poo(<?) can be any unit 
quaternion. (Recall that the matrix M is symmetric.) 

Considering the coefficient of x on each side of ( p.llf ), we obtain 

U{g)xU{g)- 1 = gxg- l t n . (2.12) 

By taking x = j, k one deduces that U (g) = p w {g)-g where p w (g) is real. Hence ADHM 
data are G-invariant if 

Poo{g)Lp w {g)~ x g~ x = L, p w (g)gMp w (g)~ 1 g~ 1 = M, (2.13) 

where p w is a real n-dimensional representation of G and is a quaternionic 1-dimensional 
representation of G. 

The reason for the notation (p w ,poo) is that A(x) can be viewed more invariantly as 
an H-linear map W — > W <8> H © Eoo where W is an n-dimensional real vector space and 
Eoo is the fibre at oo of the bundle carrying our SU(2) gauge potential. (The space W can 
in turn be identified with the zero-modes of the coupled Dirac operator, but we shall not 
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need this fact.) Here we are using the S77(2)-structure to think of Eoo as a 1-dimensional 
quaternionic vector space. 

Now the ADHM construction is natural; if we have an action of G on E (covering the 
action of G by rotations as above), such that the gauge potential is G- symmetric, then W 
and automatically become representation spaces for G. The above notation reflects 
the origin of the representations p w and p^. 

The reason why the ADHM equations become tractable after symmetry is imposed is 
simply that if W is a sum of not too many irreducible representations of G, then there will 
not be too many parameters involved in the specification of L and M satisfying (|2.13| ). Our 



next task, then, is to consider the irreducible representations of G and the construction of 
G-invariant maps between tensor products of certain of these representations. 

2.2 A little representation theory 

The problem of choosing the representations p w and p^ and of constructing the invariant 
matrices ([2.1 3|) will now be considered. Let us first compare it with the analogous problem 
of constructing symmetric Nahm data (and hence symmetric monopoles [10, [12|, |I3||). The 



single most important difference between these problems is that in the Nahm case one 
knows which representation (the analogue of p w ) is going to arise. That is because one 
knows that the Nahm data form the irreducible n-dimensional representation of SU(2) at 
the end-points, and this SU(2) really does correspond to the rotation group of M 3 . Thus 
in the cited work on symmetric monopoles, the approach was to understand explicitly how 
the standard representations of SU(2) decompose under the action of G, making use of 
the invariants corresponding to the Klein polynomials. 

In the instanton case, by contrast, we do not have such information about W. Since W 
is identifiable with a space of Dirac zero-modes, one can compute the character of W as a G- 
representation space using some equivariant index theory, but the information coming from 
this is not particularly useful. Instead we exploit the fact that the representation theory 
of finite subgroups of SU{2) can be understood very explicitly. Since every irreducible 
representation of G eventually appears in the standard representations of SU(2), the two 
approaches are closely related in principle, though this may be rather cumbersome in 
practice. 

We shall now, therefore, describe the irreducible representations of the binary icosahe- 



dral group G following John McKay's famous paper |I8|. McKay observed that a convenient 
picture of the representation theory of G is given by the extended E 8 Dynkin diagram: 

3' 



4' 



2' 



, 6 arise simply by 



Each node stands for an irreducible representation of G. Here 1, . . . , u arise simpiy uy 
restriction of the corresponding representations of SU(2) to G. The other representations 
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2', 3', 4', of dimensions 2, 3, 4 respectively, will be described shortly. First we shall explain 
the role of the edges in the diagram: the node a is joined to f3 if and only if a is a summand 
in the decomposition of f3 © 2 into irreducible representations of G. In general, this would 
lead to a directed graph, but for subgroups of SU(2), a occurs in (3® 2 iff (3 occurs in a® 2. 
Moreover, the multiplicity is always or 1. Thus we read off, for example 6© 2 = 3' ©4' ©5. 
Since moreover as S'L r (2)-representations, 6 © 2 = 5 © 7, we find that 

7 = 3' © 4' as representations of G. 

This decomposition corresponds to the existence of the icosahedral Klein polynomial of 
degree 12. Indeed the decomposition is equivalent to the existence of a non-trivial G- 
map 7 — > 7, or equally to an invariant element in 7 © 7. Now the latter contains the 
representation 13 which contains the above-mentioned Klein polynomial. 
Notice also that 4' = 2' © 2. 

In view of this last observation, it remains to describe 2' and 3'. In the character 
table of G, these two representations look exactly like 2 and 3, but with the sign of \/5 
changed. As was pointed out to the first author by J0rgen Tornehave, this extends to the 
representations, in the following sense. Identify G with a finite subset of points of H such 
that all the coordinates of these points lie in Q(a/5). In other words, for every g E G, the 
coefficients of i, j, k, (as well as the real part) are each of the form a + &a/5, where a and b 
are rationals. (We shall show how to do this in a moment.) Consider the 'conjugation' on 
Q(\/5), a I— > a' which takes a + 6a/5 to a — b\/5. Then 2' is the representation x h- > g'x, 
3' is the representation x \— > g'x(g')~ 1 (on pure imaginary x) and 4' is the representation 
x \— > g'xg^ 1 . This recipe defines a representation because of the property a'(3' = (a/3)' 
for any a, (3 E Q(v^>) arid it is clear that the character of the 'primed' representation is 
obtained by changing the sign of a/5, just as required. It so happens that only the primed 
versions of 2 and of 3 are new representations. Remark that despite the notation, 4' is not 
the primed version of 4; we hope that no confusion will result from this notation. 

In this discussion note that 3' and 4' are real representations: they act on Im(H) or H, 
viewed as 3 and 4-dimensional real vector spaces. By contrast, 2' is essentially complex; 
it is not the complexification of any real representation. We point out further that 4' has 
the following alternative description. A 5-dimensional representation of G arises through 
the identification of G with the alternating group A 5 and letting this act by permutations 
of the coordinates in R 5 . This is the sum of a trivial representation and an irreducible 
4-dimensional representation on V C R 5 , V = {y E R 5 : 2/1 + 2/2 + 2/3 + 2/4 + 2/5 = 0}. This 
representation on V is isomorphic to 4'. 

Coxeter, in 0, gives a description of the binary icosahedral group (his notation is 
(5, 3, 2)) which has the property mentioned above, namely that all coordinates lie in Q(\/5). 
For this, the icosahedron is oriented in such a way that the coordinate axes in R 3 pass 
through edge-midpoints. In particular, the half-turns about the coordinate axes lie in G. 
In terms of quaternions, these half-turns correspond to i, j, k. Then G is generated by 
these together with one other element of order two such as U2 = —(i + rj — r~ l k)/2. 
Setting 

Ui = j, U 2 = -(i + rj - r- l k)/2, U 3 = % (2.14) 
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one gets a set of generators of G. In terms of these, the five-fold rotation is given by 
A = U1U2, the three-fold rotation by B = U 2 U^ and the two-fold rotation by C = U3U1 
and 

A 5 = B 3 = C 2 = ABC 

and this equals — 1 in G but 1 when projected to G. (||, p78, eqn (7.54) and p. 69 eqn 
(6.65)). Here r = (Vo+ l)/2 is the golden ratio. 

We can now write explicitly the action of the generators U\, U 2 , U3 in the irreducible 
representations of G. We shall confine ourselves to the cases needed in this paper. Where 
necessary, we shall denote by p a the action of G in the representation corresponding to a. 
The following identities involving r will be used without comment in what follows: 



(V5-l)/2, t' = -t-\ 



T — T 



1, t + t- 1 = y/E, 



r 2 + t- 2 = 3. 



Thinking of 2 and 2' as quaternionic 1-dimensional representations, we have 

P2W = P2'(i) = hp2(j) = P2'(j) =j,p 2 (k) = p 2 '(k) = k, 

and 



P2(U 2 ) = --{% + TJ-T 1 k), P2 ,(U 2 ) 



In terms of quaternions, 3 is obtained by letting G act by conjugation on the 
imaginary quaternions Im(H). Identifying Im(H) with H 3 via the coordinates 
a\i + a 2 j + a 2 k, we find 



P3{i) = P3'W 





while 



Ps(U 2 



( 1 



— T 
-1 




-1 



P3>(U 2 



1 T 

r- 1 r 

-T 1 



-T 



.-1 



Identifying H with R 4 via the coordinates ao + a,\i + a 2 j + a^k, we obtain the action 
in the representation 4': 



Mi J 



(\ 

1 

0-1 

Vo 



\ 




-1/ 



/l \ 



<P4>(j) 



0-10 
1 

Vo 








■1/ 



and 



Pi'(U 2 



/ -1 V5 -y/h -v / 5\ 
3 1 1 

-y/E 1 -1 3 

V-v^ 1 3 -1 / 



8 



2.3 The ADHM description of an icosahedral 7-instanton 

To construct the ADHM data for an icosahedral instanton of charge n, we must choose a 
real representation of G of dimension n (the space W) and a 1-dimensional quaternionic 
representation (the space E^). Having done so, we write down the most general G-invariant 
maps and attempt to use them to solve the ADHM constraints. 

Notational remark In the rest of this paper, we shall use the term 'G-map' to denote 
any linear map between representation spaces of G that intertwines the G-actions. We 
shall also often use Schur's lemma without comment. 

Following the monopole situation, we suppose W = 7 = 3' © 4'. From McKay's 
correspondence we have 

3' © 2 = 6, 4' © 2 = 6 © 2'. 

Since L is required to be a G-map from W ®2 into E^, there is only one possibility for 
i^oo which allows L ^ 0: the only 2-dimensional representation that occurs in W © 2 is 2'. 
Hence we take 

W = 3' ©4', E^ = 2'. (2.15) 

The quaternionic matrix M is naturally viewed as a map W — > W © 2 © 2. The McKay 
correspondence can be used also to calculate the space of such G-maps. Now M breaks 
up naturally into a real component and a pure imaginary component, corresponding to 
2 © 2 = 1 © 3. The real component of M must be a multiple of the identity on the 
irreducible summands 3' and 4', yielding two free parameters. The pure imaginary part of 
M gives maps 3' — > 3' ©3, 3' — > 4' ©3, and 4' — > 4' ©3. From the McKay correspondence, 
there are no non-zero maps 3' — > 3' © 3, but there is in each case a one-dimensional space 
of maps 3' -> 4' © 3 and 4' -> 4' © 3. 

To sum up, we have the following parameter-count for G-symmetric ADHM data: one 
in the component of L that maps 4' © 2 — > 2'; two from the real part of M; and three from 
the imaginary part of M. In the following lemma we compute explicitly these invariant 
maps. Although maple has been used to assist with these computations, they are quite 
straightforward to do if one first imposes invariance under the group K C G generated by 
i, j, k. Using the above formulae, 

• The most general G-map 4' © 2 — > 2' is given by any real multiple of the row-matrix 
I = (1, h J, k). 

• The most general G-map 3' — > 4' © 3 is given by any real multiple of 



I 1 


J 


k \ 





rk 









Tl 


\ r 3 


t-H 


o J 



More specifically, this means that for each g e G, we have 

g'l pi'igy^g' 1 = I, 
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and 

P4>(g)gB py{g)~ 1 g~ 1 = B. 
Taking the quaternionic conjugate of this, 

p 3/ (g)gB^p 4 ,(g)- 1 g- 1 = Bl 

It follows that 



/3 -i 




-k\ 


i 3 


k 




j —k 


3 


% 


\k J 


—i 


3 / 



satisfies p 4 i{g)g BB^ p 4 i{g)~ 1 g~ 1 = BB\ so its imaginary part gives the unique G-map 
4' <S> 3 — > 4'. Since this is skew symmetric, it cannot be used as a diagonal block in the 
matrix M from the ADHM data, since M is required to be symmetric. It follows that the 
most general (5-invariant ADHM data, with W = 3' © 4', is given by 



M 




(2.16) 



where a, b, c, d are real numbers and we have used the fact that B is pure imaginary to 
write B l = —B^. One computes 



B^B = 41 3 , lU = 41 4 - BB* 



so that 



(4a 2 + b 2 ) 1 4 + (c 2 - a 2 )BB^ c(b - d)B 
c{b-d)B^ (4c 2 + d 2 )l 3 , 

Thus the ADHM constrainsts are satisfied iff a 2 = c 2 and c(b — d) — 0. If c = 0, then a = 
and the top row L of M is identically zero. This yields a singular instanton and hence is 
not allowable. Hence c ^ 0, 6 = d and a = ±c. Since, moreover, 



(I 







u 





\ 


-1,3/ 




/ aZ \ 



the two choices of sign lead to gauge-equivalent ADHM data and hence to gauge-equivalent 
instantons. 

What we have found is a charge-7 icosahedral instanton that is unique up to the ob- 
vious freedom to translate along the x 4 axis (the parameter b) and overall scale ((a, b) \— > 
(Aa, \b)). If we centre the instanton at the origin of M 4 then the ADHM data are given 
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by any real multiple of 



/l 


i 


J 


k 

























i 


J 


k 



















rk 


















T- l k 





T% 
















rj 


t-H 







i 





T- l k 


rj 













J 


rk 





t-H 













\k 




Tl 














/ 



(2.17) 



A maple calculation shows that A(x)^A(x) is invertible for every x, so that this is a non- 
singular instant on. 



2.4 The associated Skyrmion 

Having constructed the ADHM data for an icosahedral instanton we now wish to make 
use of this to compute a Skyrme field. Recall the proposal of Atiyah and Manton[0 which 
results in the following explicit prescription for the Skyrme field 

U(x) =Vexp (- J +0 ° A 4 (x,x 4 ) dx^j. (2.18) 

Here £/(x) is the SU (2)- valued Skyrme field in K 3 , V denotes path ordering and is the 
gauge potential of a Yang-Mills instanton field in R 4 . 

For an instanton of charge n this holonomy produces a Skyrme field with baryon number 
n. Although this procedure does not give exact solutions to the Skyrme model it does give 
fields which are good approximations to important Skyrmion configurations, in the sense 
of having not only the correct symmetries but also energies which are only a few percent 
above those of the numerically known Skyrmion solutions. For example, the 1-instanton 
generates a hedgehog Skyrme field and by adjusting the scale of the instanton, which 
may be regarded as a free parameter in the approximation, it is possible to obtain an 
approximation whose energy is only 1% above that of the true solution (which is known 
only numerically). 

The minimal energy Skyrmions of charge two, three and four have axial, tetrahedral 
and cubic symmetry respectively]?]]. In each of these cases instantons have been found 
with the correct symmetries, so that computing their holonomies produces Skyrme fields 



which are good approximations to these Skyrmions |2|, [TS|. It should be noted that in these 
multi-instanton examples the holonomy can not be computed analytically and therefore 
numerical methods must be employed. Specifically, if U (x, x 4 ) denotes the solution of the 
matrix ordinary differential equation 

dJj = -Ajj (2.19) 

with the initial condition £/(x, — oo) = 1 2) then U(x) = U (x, oo). The above set of ordinary 
differential equations, with x regarded as a parameter, need to be solved numerically, for 
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which we employ a standard Runge-Kutta method. In principle, since we have the ADHM 
data explicitly, it is possible to obtain an exact analytic expression for A 4 by performing 
some quaternionic linear algebra. However, this is a non-trivial computation, even with 
the use of a symbolic computer algebra package, and is of little practical use. Furthermore, 
since we are employing a numerical algorithm to solve equation ( |2.19| ) then it is sensible 
to compute A4 numerically also. To achieve this, note that equation (|2.4]) states that the 
vector N(x) is orthogonal to each of the columns of the matrix A(x). Thus we can compute 
N(x) using a quaternionic Gram-Schmidt orthonormalization process and hence A&(x) by 
using a finite difference approximation to equation (|2.5| ). 

Applying our numerical scheme to the ADHM data of the icosahedral 7-instanton de- 
rived above we obtain a Skyrme field whose baryon density isosurface is displayed in Figure 
1.11. This surface, which looks identical to that obtained from full field simulations of the 
Skyrme model 0], resembles a dodecahedron, with the baryon density being maximal at the 
vertices of the dodecahedron and holes at the centre of each face. We have not attempted 
to find the instanton scale at which the energy of this Skyrme field is minimized since to 
perform an accurate calculation would require a more substantial amount of computing 
time. The important point is that we have demonstrated that an icosahedrally symmetric 
7-instanton exists whose holonomy produces a good approximation to the minimal energy 
7-Skyrmion. 



3 Seven Skyrmion scattering 

The aim of this section is to obtain a family of instantons which describes a seven 
Skyrmion scattering process in which seven well separated Skyrmions merge to form the 
dodecahedral 7-Skyrmion. In order to understand this scattering it is useful to first consider 
the analogous situation for monopoles, since from numerical simulations it appears that 
many Skyrmion scattering events are remarkably similar to known monopole scatterings |J. 

A dodecahedral 7- monopole exists |12| whose energy density isosurface looks very similar 
to the baryon density isosurface shown in Figure 1.11. In fact it was known earlier that 
a tetrahedral 3-monopole and a cubic 4-monopole also exist [|H]] and again they resemble 
the corresponding Skyrmions fl3f] . In studying symmetric monopoles it is useful to consider 
rational maps, which we outline below. 

The n-monopole moduli space is diffeomorphic to the space of degree n rational maps 
between Riemann spheres, with the equivalence relation that two maps that are equal after 
a rotation of the target sphere are identified. The rational map arises from the monopole 
as the scattering data along a half-line emanating from a chosen origin[|TJ[]. Explicitly, let 
z be a point on the Riemann sphere and consider Hitchin's equation 

(A-i$)s = (3.1) 

for the two-component field s, along the radial half-line through the point z. Here D r is 
the covariant derivative in the radial direction and $ is the Higgs field. Up to a constant 
multiple, there is a unique solution s = (si,S2)* which decays as r — > 00. Let R be the 
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ratio of the components of this solution evaluated at the origin, that is, R = j*\r=o- Now 
consider how R varies as we choose a new direction for the half-line by changing the value 
of z. Then, as proved by Jarvis[14], R is a holomorphic function of z of degree n, where n 



is the charge of the monopole fields occuring in (|3.1| ). The effect of a gauge transformation 
is to transform R by an £77(2) Mobius transformation, and after taking this equivalence 
into account, there is a one-to-one correspondence between n-monopoles and rational maps 
R(z) of degree n. 

If a map R(z), of degree n, is G-invariant (up to Mobius transformations) then there 



is an n-monopole with symmetry G, and vice versa. In ref.[]Il] many symmetric maps are 
presented but the one of relevance here is the following degree seven map which arises after 
the imposition of the symmetry Th 

bz 6 - 7z 4 - bz 2 - 1 

R ( z ) = -r-« — r~A — =-5 — 7T- 3.2 

v ; z(z e + bz 4 + 7z 2 - b) v 1 

Th is the group of rotations of a tetrahedron extended by inversion symmetry and, after a 
choice of orientation, the above one-parameter family, with b real, gives all such maps. Since 
this one-parameter family is the fixed point set of a group action then it is a geodesic in the 



7-monopole moduli space. Using the geodesic approximation [|I7| this family describes a 
low energy seven monopole scattering process as b varies along the real line from — oo to oo. 
Changing the sign of b can be undone with a Mobius transformation plus the replacement 
z i — > iz, which corresponds to a rotation by 90° about the xi-axis. lib = ±7/y/E then 
the map has icosahedral symmetry, and represents the dodecahedral 7-monopole and its 
dual, whereas at b = 0, which is the midpoint of the scattering process, the map has cubic 
symmetry. In the limit as b — > oo the map degenerates to R{z) = z, which represents a 
single monopole at the origin, the other six monopoles having moved off to infinity along 
the Cartesian axes. 

In summary this geodesic models a scattering event where six monopoles, moving in 
along the Cartesian axes, merge with a single monopole at the origin to form first a 
dodecahedron and then a cube, after which the process reverses but with a 90° rotation. 
The purpose of the remainder of this section is to present a one-parameter family of Skyrme 
fields which describe a similar scattering of Skyrmions. At this point it is important to 



note that although an ansatz exists for Skyrme fields in terms of rational maps|Tl|, this 
approximation only works for Skyrme fields which have a shell-like structure. Thus, for 
example, using the rational map ( |3.2| ) with 6=7/ \/5 gives a good approximation to the 
dodecahedral 7-Skyrmion, but the rational map ansatz breaks down for large b and can 
not be used to describe well separated Skyrmions. Thus we need to turn to the instanton 
approximation to attempt to produce Skyrme fields which describe this process. 

The upshot of the above discussion is that we now want to consider 7-instantons with 
symmetry Th- 
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3.1 Tetrahedral deformations of the icosahedral instanton 



We have described the general framework for the construction of symmetric ADHM data. 
What we seek now is the most general family of tetrahedral 7-instantons that contains the 
icosahedral 7-instanton. Denote by T the binary tetrahedral group, i.e. the double cover 
in SU (2) of the rotation group of the tetrahedron. To study tetrahedral deformations, we 
realize T as a subgroup of G generated by i, j, k and the three-fold rotation (1 + i + j + 
k)/2. (We continue to take G to be generated by the C/j of (|2.14Q .) Conjugation by this 
quaternion gives a 120°-rotation about the axis in the direction (1, 1, 1) in H 3 . Then the 
given representations of G yield representations of T and we can attempt to follow the 
procedure described above. It is preferable first to impose the inversion from the subgroup 
Th, however, since this forces the two diagonal blocks in M to be zero. 



To explain this, we must first show in what sense the data (|2.17f) are inversion sym- 
metric. For this we must again find compensating gauge transformations J\ and J 2 such 
that 

A(-x) = JiA(x)J 2 



r 2 \ 



It is easy to see that the essentially unique choice for this, in the case of the ADHM data 
( p. 171 ), is given by 



f-1 


V 





14 




J-2 



-14 







la 




The reason for the simple form of these matrices is that the inversion x 1— > — x is central 
and so must act by ± the identity in any irreducible representation of G. 

Now any ADHM data that is inversion-symmetric in this sense must have the form 



M 



relative to the usual block-decomposition. The row-vector v and 4x3 matrix C are 
constrained by having to be symmetric under T, where this acts by restriction of the 
given G-representations. Since 2 = 2' and 3 = 3' as T representations, we have that 
4' = 2(g>2 = l©3as T-representations. One may check that v is invariant iff it has the 
form (Aq, Aii, Xij, X±k) and that C is invariant if it has the form 



C 



where A , Ai, A 2 , p, q, r are real numbers. 

The next task is to impose the ADHM constraints; these yield the equations 



/A 2 2 


A 2 j 


X 2 k\ 


p 


qk 


rj 


rk 


P 


qi 


V q j 


ri 


P / 



lm{v^v + CC l ) = 0, Im^C) = 0. 
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It is straightforward to show that these reduce to the three equations 

A Ai = X 2 (p + q — r), X\ = qr - p(q - r), X\ = qr + p(q - r). 

Solving for the A's in terms of p, q, r, we conclude that the general Th invariant ADHM 
data are given by 



M 



f A 


Xii 


Au 


Xik 








\ 














X 2 i 


hj 


X 2 k 














P 


qk 


rj 














rk 


P 


qi 














QJ 


ri 


P 


X 2 i 


p 


rk 


QJ 











A 2 j 


qk 


P 


ri 











\X 2 k 


rj 


qi 


P 








o ) 



(3.3) 



where 



Ao = ip + Q - r) 



qr + p(q 



\ qr — p(q 



Ai 



qr — p(q — r), A 2 



qr + p(g 



(3.4) 



Here we have three free parameters, p,q,r, but we require only a one-parameter family 
for our application. One of the free parameters is an overall scale factor, which corresponds 
to the multiplication of M by a scale, and this can only be determined by computing the 
energy of the associated Skyrme field. We will therefore set this scale to one, as it can be 
easily reintroduced later. As p, q, r are homogeneous coordinates we can fix this scale by 
setting r = q~ x . 

To get some insight into the interpretation of the parameters, let us consider the special 
case p = 0. Then from ( |3.4D , Ao = (q — q^ 1 ), X\ = X 2 = 1; denote the corresponding matrix 
( p.3| ) by M(q). If q — r, then Ao = 1 and M(r) coincides with the icosahedral matrix 
( p.!7| ). If q is interpreted as the exponential of a 'time parameter' t, then q runs from 
to oo as t runs from — oo to oo, and time-reversal is the transformation q — > This 
transformation is equivalent to a 90° rotation about the X\ axis. Indeed if we denote by 



N(q) the ADHM data obtained from M(q) by replacing q by q 
—j, then we have 



N(q) 




M{q) 



Ri 






R2 



by i, j by k and k by 



(3.5) 



if 



Ri 



f 1 




Vo 





-1 










1 






-1 














, R2 = 


(i 











-1 




J 





(3.6) 



We conclude that the value q 
by a 90° rotation) and that at q 



T 



~ l is also a dodecahedral instanton (obtained from ( |2.17|) 
we have an instanton symmetric under the symmetry 
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group of the cube (since this group is generated by the tetrahedral group together with 
any 90° rotation about one of the coordinate axes). 

This one-parameter family thus has all the symmetry properties expected of the seven 
Skyrmion scattering process. However it turns out that it does not have the correct asymp- 
totic behaviour as q goes to or oo. This was verified both by numerical work and asymp- 
totic analysis. The problem can be traced to the fact that the first entry of the top row 
(q — q^ 1 ) blows up as q — > 0, oo. 

Therefore we try to find a more general family with the same symmetry properties but 
allow p 7^ to improve the asymptotic behaviour. In particular we want time-reversal to 
continue to correspond to the replacement q h- > q~ x . Then it can be seen that if we make 
the replacements 

q i— > q^ 1 ,p I— > —p,i i— > i, j i— > k, k i— > —j (3-7) 
which, because of ( |3.4j) , result in 

A >->• -A , Ai h-> Ai, A 2 i-> A 2 , (3.8) 

equivalent ADHM data are obtained. (The 'compensating gauge transformation' is as in 
( p.5| ) and ( |3.6Q .) Note that the fixed point set of the time reversal transformation, in other 
words the midpoint of the scattering process, is given by q — 1, p — 0, which indeed has 
cubic symmetry as it should. 

These arguments show that if p is any function of q with the properties 

p( g -i) = -p(q), p( r ) = (3.9) 

then the corresponding family of ADHM matrices will have the correct symmetry prop- 
erties: the first condition gives that the data at q and at q~ l are equivalent up to a 
90°-rotation; the second ensures that the data at r reduce to the icosahedral data fl2.17| ). 

We are now left with finding the variable p as a function of q; this can not be determined 
from symmetry arguments alone. A simple function satisfying (|3.9|) is given by 



, = _(«-«-')«« -f)'-!). (3 . 10 ) 

(,_, -1)4 + 1 ^ 

The numerator of p is the simplest function having the required zeros and symmetry 
properties, whereas at this stage the only fact we know about the denominator is that it 
must be a symmetric function of (q — q~ x ). However, as we shall now see, the form of 
the denominator is highly constrained by examining the asymptotic limit of the scattering 
process, where all seven instantons are well separated, and this leads naturally to the given 
solution. 

The asymptotic out state corresponds to the limit q — ■> oo, and in this limit the leading 
order behaviour of p is 

p = -q- 1 + 0(q- 5 ). (3.11) 

Hence in this limit we find that 

A = ^= + 0(q- 2 ), X, = V2 + 0(q- 2 ), A 2 = q- 1 + 0(q~ 3 ). (3.12) 
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Neglecting negative powers of q we thus arrive at the asymptotic ADHM data 



(1/V2 


y/2i 


y/2j 


V2k 








o\ 






































qk 























qi 














QJ 

















qj 














qk 

















V 





qi 











0/ 



(3.13) 



After a gauge transformation by the matrix 

/V2 



T 



1 

V2 












1 





1 







1 





1 







1 
1 










1 

-1 








1 










1 







(3.14) 



■1/ 



we obtain 




T 



(I 


i 


J 


k 


k 


i 


j\ 


























—qk 























—qi 























-qj 























qj 























qk 





Vo 

















qi) 



(3.15) 



This matrix is of the form identified by Christ, Stanton and Weinberg!! as representing 
well separated instantons. The leading order terms lie on the diagonal of the square part 
of the matrix and determine the instanton positions, giving one instanton at the origin and 
the other six on the Cartesian axes at a distance q from the origin. The terms of next order 
all lie on the top row of the matrix and these give the scales and SU(2) orientations of the 
instantons. The fact that all the entries on the top row have unit length means that each 
instanton has scale one. Since the first entry in the top row is one, then on computing the 
holonomy the Skyrmion at the origin will be in standard orientation, whereas, for example, 
the fact that the second and sixth entries in the top row are i, means that the Skyrmions 
located on the a^-axis have an orientation which is obtained from the standard one by 
rotating the Skyrmion by 180° around the xi-axis. These kinds of configurations, that is 
three collinear Skyrmions such that the outer Skyrmions are rotated by 180° about a line 
perpendicular to the line joining the inner Skyrmion, are known to give attractive initial 
conditions for full field simulations ||. Thus we see that our ADHM data has the correct 
asymptotic properties to produce the configuration which we require. 
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Now we address the possible freedom in the choice ( |3.10| ) that we have made. The 
main effect of changing the denominator in ( p,10|) is to alter the scale of the instanton at 
the origin in the limit in which the other six are far from it. The requirement that this 
scale is finite as q —>■ oo determines that the leading term in the denominator of p must be 
(<? — Q ) w ith coefficient one. Only even powers of (q — q^ 1 ) are allowed by symmetry and 
it can be shown that the coefficient of (q — q^ 1 ) 2 must be zero if the scale of the instanton 
at the origin is to be not only finite but equal in value to the scale of the instantons which 
are moving along the Cartesian axes. There remains the freedom to change the value of 
the constant in the denominator of ( [3.10|) but clearly this has little effect since it is only 
relevant for small values of (q — q^ 1 ) and the numerator contains (q — q^ 1 ) as a factor. 

Using this ADHM data for increasing values of q we compute Skyrme fields whose 
baryon density isosurfaces are shown in Figure 1 . The various values of q corresponding to 
each figure are given in the table below. 



Fig 


1 


2 


3 


4 


5 


6 


7 


8 


q 


0.40 


0.44 


0.48 


0.50 


0.62 


0.77 


0.87 


1.00 


Fig 


9 


10 


11 


12 


13 


14 


15 




q 


1.15 


1.30 


1.62 


2.00 


2.10 


2.25 


2.50 





Table 1. Parameter values for the scattering shown in Figure 1. 

From Figure 1 we see that indeed our family of instantons describes the sought after 
scattering process. Figure 1.1 clearly shows the six Skyrmions on the Cartesian axes and 
a Skyrmion at the origin. As the Skyrmions approach, Figure 1.2, the one at the origin 
shrinks until it disappears completely, Figure 1.3. The Skyrmions then merge until the 
dodecahedron is formed, Figure 1.5, after which the configuration twists until it turns into 
a cube, Figure 1.8. This process is then reversed but accompanied by a 90° rotation around 
the xi-axis, so that for example the dual dodecahedron is formed, Figure 1.11, and the 
Skyrmions finally separate again along the Cartesian axes, Figure 1.15. 

An important difference between monopoles and Skyrmions is that monopoles are BPS 
solitons, and therefore all configurations of the same charge have equal energy, whereas for 
Skyrmions this is certainly not the case and the potential energy of different configurations 
is an important factor in considering the dynamics. We have not computed the scale factor 
of each configuration to minimize the energy of our Skyrme fields since, although this could 
be done if required, it would involve a substantial amount of computing time. However, we 
have ensured that the scale factor is the order of unity for all q, so that we get an accurate 
representation of the baryon density isosurfaces. Qualitatively we know that for all the 
family of Skyrme fields the energy per Skyrmion is less than it is for seven well separated 
Skyrmions, and that the minimum energy configurations are the two dodecahedra. From 
numerical simulations it is known that the energy of the cubic 7-Skyrmion is above 
that of the dodecahedron but it is still substantially less than that of seven well separated 
Skyrmions. Thus the true dynamical evolution depends upon the initial speeds of the 
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incoming Skyrmions, together with the amount of energy lost through radiation as the 
process evolves. However, it is expected that if the incoming speeds are great enough then 
the whole scattering process displayed in Figure 1 would take place. Radiation effects will 
mean that for most speeds the incoming Skyrmions will eventually get trapped at one 
of the dodecahedra, and perhaps if the Skyrmions are initially static then only the first 
portion of the scattering process will occur and the cube may never be formed, but this 
depends upon the amount of radiation generated. By performing full field simulations, 
using the numerical code described in || ||, with initial conditions given by the instanton 
generated Skyrme field, it has been verified that the true dynamical evolution does follow 
the sequence described above and the family of instanton generated Skyrme fields provides 
an accurate approximation to the Skyrmion scattering process. 

In current approaches to the quantization of Skyrmions a first step is to examine the 
vibrational modes of the minimal energy Skyrmion at each charge [Q. Thus for charge seven 
it is the vibrational modes of the dodecahedral Skyrmion which need to be studied and 
clearly one mode is the tetrahedral deformation we have displayed. Thus if the amplitude 
of the deformation is sufficient then one of the important vibrational modes will be the 
one considered here where the dodecahedron deforms to its dual via a cube. 



4 Conclusion 

We have used the ADHM construction to present a charge seven instanton with icosahe- 
dral symmetry whose holonomy generates a Skyrme field which approximates the minimial 
energy dodecahedral 7-Skyrmion. Furthermore, by imposing tetrahedral symmetry we have 
found a family of ADHM data which we used to generate Skyrme fields that model a seven 
Skyrmion scattering process that results in the formation of the dodecahedral 7-Skyrmion. 

There are a number of other highly symmetric Skyrmions, such as the icosahedrally 
symmetric charge seventeen Skyrmion which resembles a buckyball ||ll|| , and the methods 



used here could also be applied to construct the ADHM data of the predicted correspond- 
ing instantons. In particular it would be interesting if the instanton approach led to an 
understanding of the formation of the buckyball 17-Skyrmion from individual Skyrmions. 
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Figure Caption 

Fig. 1. Baryon density isosurfaces for a family of charge seven Skyrme fields obtained from 
7-instantons with tetrahedral symmetry. 
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